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ABSTRACT 

We formulate the first order Fermi acceleration in shock waves in terms of the random 
walk theory. The formulation is applicable to any value of the shock speed and the 
particle speed, in particular, to the acceleration in relativistic shocks and to the in- 
jection problem, where the particle speed is comparable to the fluid speed, as long as 
large angle scattering is suitable for the scattering process of particles. We first show 
that the trajectory of a particle suffering from large angle scattering can be treated 
as a random walk in a moving medium with an absorbing boundary (e.g., the shock 
front). We derive an integral equation to determine the density of scattering points of 
the random walk, and by solving it approximately we obtain approximate solutions 
of the probability density of pitch angle at and the return probability after the shock 
crossing in analytic form. These approximate solutions include corrections of several 
non-diffusive effects to the conventional diffusion approximation and we show that 
they agree well with the Monte Carlo results for isotropic scattering model for any 
shock speed and particle speed. When we neglect effects of 'a few step return', we 
obtain 'the multi-step approximation' which includes only the effect of finite mean 
free path and which is equivalent to 'the relativistic diffusion approximation' used by 
Peacock (1981) if the correct diffusion length is used in his expression. We find that 
the multi-step approximation is not appropriate to describe the probability densities 
of individual particles for relativistic shocks, but that the pitch angle distribution at 
the shock front in steady state is in practice quite well approximated by that given 
by the multi-step approximation because the effects of finite mean free path and a 
few step return compensate each other when averaged over pitch angle distribution. 
Finally, we give an analytical expression of the spectral index of accelerated particles 
in parallel shocks valid for arbitrary shock speed using this approximation. 

Key words: acceleration of particles - methods:analytical - shock waves - cosmic 
rays. 



1 INTRODUCTION 

It is now widely accepted that the first order Fermi acceler- 
ation, which is driven in various astrophysical shocks, pro- 
vides a factory of non-thermal energetic particles. For exam- 
ple, recent observations discovered synchrotron X-rays from 
energetic electrons with energ ies above 10 14 eV in the shell of 
the supernova remnant 1006 (Koyama et al. 1995). Since the 



basic th eory was proposed more than 20 years ago by Bell 
([L978a|), Krymski (1977), Axford et al. (1977), and Bland- 



ford and Ostriker (1978), many papers have dealt with fun- 
damental mechanisms of shock acceleration. The mechanism 
has been applied to many situations ranging from interplan- 
etary shocks to relativistic shocks in active galactic nuclei 
or gamma-ray bursts (see for review Blandford & Eichler 



1987 and Kirk & Duffy 1999). Although most of the recent 
theoretical interests have been concentrated on non-linear 
problems where the shock structure is modified by acceler- 
ated particles (Ellison, Baring & Jones 1996, Berezhko & 
Ellison 1999) , some of the linear problems in the test parti- 
cle approximation still remain to be clarified and need more 
close examinations, especially for relativistic shocks, oblique 
shocks and injection of seed particles. 

The first order Fermi acceleration in shock waves with 
test particle approximation is the most basic problem of the 
shock acceleration. As a particle gains energy whenever the 
particle repeats shock crossing and recrossing between the 
upstream and downstream of the shock front, it is essen- 
tial to calculate the average energy gain per this cycle and 
the number distribution of the repeated cycles before the 
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particles escape to far downstream. These are determined 
by the shock speed, the compression ratio of the shock, the 
pitch angle distribution of particles at the shock crossing 
and the return probability of particles, i.e., the probability 
of recrossing the shock front for particles that cross the front 
from the upstream to the downstream. As is well known, the 
energy spectrum of particles accelerated by this mechanism 
takes a power law form, F(E)dE oc E~° dE. For the non- 
relativistic shocks, the conventional diffusion approximation 
can be used to obtain the universal power law spectrum with 
an index of a — (r+2) /(r— 1), where r is the compression ra- 
tio of the shock (Bell 1978a) . However, for relativistic shocks, 
oblique shocks and acceleration of low energy particles, in 
which the fluid speed is not negligible compared with the 
particle velocity, anisotropy of particle distribution at the 
shock front becomes large and the diffusion approximation 
is not relied on. Thus, one should develop more sophisticated 
treatments beyond the diffusion approximation even in the 
linear regime. 

There are three ways to do this. One is to solve the 
transport equation with appropriate collision operators, as 
was developed by Kirk & Schneider (198 7, 19 88) for rela- 
tivistic shocks, and by Malkov & Volk ( 1995 ) for the in- 
jection problem. They solved for eigenvalues and eigenfunc- 
tions of the transport equation both for upstream and down- 
stream and by using the matching condition at the shock 
front they obtained the spectral index a and the distribu- 
tion function. This method works well when the pitch angle 
diffusion is a main process of the scattering, while it seems 
to have met some difficulties when the simpler case of pure 
large angle scattering is considered. 

Another approach is the sin gle-par ticle approach, which 
was originally used by Bell (1978a) for n on-relativistic 



shocks and was extended by Peacock ( 1981 ) for relativis- 
tic shocks. Although this method is in principle equivalent 
to the former method, it is more intuitive and more tractable 
as will be described in this paper. To det ermine the pitch an- 
gle dis tribution at the shock front, Bell ( 1978a ) and Peacock 
(1981) basically utilized the distribution functions obtained 
by ma crosc opic methods such as the diffusion equation. Pea- 
cock (1981) also used an upstream distribution function, 
which is realized far upstream when particles suffer large 
angle scattering; these functions may be different from real 
ones at the shock front as was shown in Kirk & Schneider 
(1987, 1988). 

The third one is the Monte Carlo simulation in which 
the motion of each particle is traced faithfully (Ellison, Jones 
& Reynolds 1990; Ostrowski 1991; Bednarz & Ostrowski 
1998). While Monte Carlo simulation can make a direct es- 
timation of the distribution function and the spectral index, 
a large scale simulation is needed to obtain sufficiently ac- 
curate results. Physical interpretations of simulation results 
in terms of analytical models are desirable. 

The results of acceleration in such situations in fact de- 
pend on the nature of scattering mechanisms of particles 
by magnetic irregularity existing in the background plasma. 
This is because the anisotropy of particle distribution at the 
shock front is strongly subject to the detail of the scatter- 
ing process. Especially, two conventional scattering models, 
large angle scattering and pitch angle diffusion, lead to sig- 
nificantly different spectral index in relativistic or highly 
oblique shocks (Kirk & Schneider 1988; Ellison, Jones & 



Reynolds 1990; Naito & Takahara 1995; Ellison, Baring & 
Jones 1996). It also may be the case for acceleration of low 
energy particles. The pitch angle diffusion is derived by the 
quasi-linear theory of plasma turbulence (e.g., Melrose 1980) 
but this description is adequate only for weak turbulence 
while it is often supposed that the turbulence can be strong 
in astrophysical shock environments (see Ellison, Jones & 
Reynolds 1990 and references therein). On the other hand, 
the large angle scattering, which isotropizes the particle di- 
rection in a single scattering independently of the direction 
before scattering, may mimic some effect of strong turbu- 
lence, though there is no theoretical justification for this. 
Theory of particle transport in strongly turbulent field re- 
quired for realistic investigation of this problem is very dif- 
ficult and not yet established, to our knowledge. 

In this paper, we adopt the large angle scattering model 
regarding that it mimics the scattering under strongly tur- 
bulent field, and investigate the nature of the shock accel- 
eration in highly anisotropic situation, basically taking the 
second approach mentioned above. The specific large angle 
scattering model used in the paper is prescribed in Section2, 
and assumes that the particle direction after scattering does 
not depend on the direction before scattering (although this 
feature is not always true for general large angle scatter- 
ing, we use the term 'large angle scattering' for the specific 
model in this paper). This mathematical simplicity allows 
us to treat the problem of particle transport as a random 
walk. Thus, our approach is based on the theory of random 
walk in a moving medium not on the diffusion theory to ex- 
amine the pitch angle distribution at the shock front and the 
return probability of the particles from the downstream of 
the shock. Therefore, it can be applied to any value of the 
shock velocity. The paper is organized as follows. In Sec- 
tion ^| we show that the motion of a particle suffering from 
the large angle scattering can be treated as a random walk 
in a moving medium, and formulate the random walk using 
the probability theory. In Section we derive analytical ap- 
proximate solutions, and devote to a specific case of isotropic 
scattering in Section ^. Finally, we apply these results to the 
shock acceleration in Section pi 



2 METHOD 

In this paper, we treat the first order Fermi acceleration in 
parallel shocks for any value of the shock speed. We assume 
the test particle approximation and adopt large angle scat- 
tering as the scattering process of particles. 

In order to know how particles are accelerated, the prop- 
erties of particle trajectories are to be investigated in detail. 
When particles move around on each side of the shock front, 
their motion can be regarded as one-dimensional because, in 
plane-parallel shocks, only the position along the field line 
is relevant. If the scattering process is assumed to be de- 
scribed by large angle scattering, the motion of particles on 
each side of the shock front can be treated as a random walk 
in a moving medium as is described in this section. Thus, 
the problem will be reduced to the random walk of particles 
in a moving medium with an absorbing boundary; the shock 
front becomes such a boundary both for upstream and for 
downstream. 

In this section, first, we formulate a random walk treat- 
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ment of the motion of particles following the large angle 
scattering and examine the properties of the random walk 
in a moving medium. Although the method described in this 
section may also be used for general problems of particle mo- 
tion in a moving medium if this scattering model is suitable, 
we confine our attention to the Fermi acceleration in shock 
waves in this paper. 

In the following argument, we assume that the velocity 
of scattering centres is uniform on each side of the shock 
front and is be equal to the velocity of fluid (therefore, we 
call the velocity of the scattering centres the fluid velocity 
in the following). It is also assumed that energy loss mech- 
anisms of particles (e.g. synchrotron loss) is negligible, and 
that the scattering of particles is elastic in the fluid frame, 
i.e., the particle energy does not change upon scattering. 

Hereafter, we take the unit c = 1 (c is the speed of 
light) and take the spatial coordinate x along the magnetic 
field line or the shock normal toward downstream direction 
so that the fluid velocity, V{, is positive on each side of the 
shock front. We use subscript (+) and (-) to denote the 
downstream direction (+a;) and upstream direction (—a;), 
respectively. Physical quantities measured in the fluid frame 
are expressed with a prime. 

2.1 Large angle scattering model 

We prescribe the large angle scattering model in the fluid 
frame by the following three conditions. 

(i) The energy of particle measured in the fluid frame is 
conserved at scattering. Thus, while particles stay in either 
the upstream or the downstream region, particle speed v' 
does not change in the fluid frame. 

(ii) The probability density of displacement of the particle 
along the magnetic field line between successive scatterings 
is assumed to obey an exponential distribution with mean 
free path A = X(v' , //). (Here, the term 'the mean free path' 
means the displacement along a magnetic field line and not a 
length on the orbit of its gyro motion.) Then, the probability 
density of the displacement Ax' is written as 

p(Ax';v',fj,')dAx' = ( j^^dAx' (// • Ax' > 0) (1) 
[ (otherwise), 

where the signs of Ax' and \x must be the same. 

(iii) The pitch angle cosine of particle [s! after scattering 
is determined according to the probability density function 
P fl i(fj,';v') independently of the pitch angle before scatter- 
ing. Below, for simplicity, we represent this function as P M / . 
Clearly, it must satisfy J P^d/j, = 1. 

The functional form of X(v',fi') and P M < (// ;v') still is free 
and various scattering processes can be adopted. An often 
used case is that the mean free time is independent of [i! and 
the scattering is isotropic, which we consider in Section ^. 

Because the problems of particle transport with a 
boundary are to be considered in the following, it is bet- 
ter to describe the position of the particle in the reference 
frame in which the coordinate of the boundary position takes 
a fixed value, and scattering centres are flowing with fluid 
speed Vf. We call this reference frame the boundary rest 
frame. (For the shock acceleration, this frame is of course 
the shock rest frame for each side of the shock front. But we 



consider general case here.) Therefore, below, we describe 
the speed v' and pitch angle cosine fjf of particles in the fluid 
frame and describe the position of particle x and the time t 
in the boundary rest frame. (Since we consider only steady 
states below, we need not consider about t in the following 
in fact.) This treatment is widely used for various transport 
equations such as diffusion convection equation or radiation 
hydrodynamics (Kirk, Schlickeiser & Schneider 1988). Thus, 
we describe the scattering law in terms of the displacement 
measured in the boundary rest frame, Ax, which is given 
through a Lorentz transformation of Ax' as 



Ax = r f (Ax' + VtAt') =r f Ax'(l + — ) 

v' ' \l' 

= r f Aa/(l + ii) 

where v is defined as 
Vt 



(2) 



(3) 



and Ff = 1 — Vf 2 . Here, v represents the degree of ad- 
vection effects to the particle motion in the boundary rest 
frame and is a fundamental parameter in the random walks 
in a moving medium described below. Clearly, if v > 1, a 
particle cannot advance against the flow of scattering cen- 
tres, which is not of interest for the shock acceleration. 
Hence, we consider only the case of v < 1 below. 

It can be shown easily that when a particle moves to- 
wards (+) direction Aa; > 0(/i > 0) in the boundary rest 
frame, the range of the pitch angle cosine in the fluid frame 
is — v < fj,' < 1, and for (-) direction Ax < 0(/i < 0), the 
range becomes —1 < fjf < —v. Then, the probability distri- 
bution of Ax becomes 



p(Ax; v , n)dAx = 



~dAx 



where 

A* = A>V) := r f A(t/,/i')|l+-^l- 

A 1 



((// + v) ■ Ax > 0) (4) 
(otherwise) 



(5) 



Here, A* denotes the mean free path of the particle measured 
in the boundary rest frame and can depend on [i even if v' 
is the same. For later argument, we define the following two 
scale lengths, 



L+ := max(A*; — v < p! < 1), 
L_ := max(A*; — 1 < fjf < —v) 



(6) 



where L+ and L_ denote the scale of scattering length in the 
boundary rest frame for (+) direction and for (-) direction, 
respectively. 

It should be noted that, even if both A and P^i depend 
on v' , when the dependences are separable into v' and //, 
the properties of particle motion is determined only by v 
except that the temporal and spatial scales can be different. 
As far as only the probability of return and the pitch angle 
distribution at the shock front are to be considered, these 
scales are not relevant at all. Then, one can choose the scale 
of the time and space in order for the computation to become 
the simplest. 
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2.2 Random walk of particles in a moving medium 

In the large angle scattering model described in the previous 
subsection, since the pitch angle of a particle after scattering 
is determined independently of that of before scattering and 
v is conserved, the displacement of the particle between suc- 
cessive scatterings is determined independently and equally 
on each scattering. Because of this property, the particle 
motion can be treated as a well-defined random walk. In 
this subsection, we formulate such a random walk treatment. 
Our arguments as to the random walk in this subsection are 
based on the Chapter 2 of Cox & Miller (1965). 



2. 2. 1 Random walk for large angle scattering models 

In the random walk considered here, the displacement of 
the particle Ax takes a continuous value rather than a dis- 
crete value. For such a random walk, one of the fundamen- 
tal concepts is the probability density function (p.d.f.) of 
the displacement for each step of random walk, f(Ax). For 
the large angle scattering model, we can write this function 
in terms of p(Ax; v' , //) and P^t, which introduced in the 
previous subsection. Using equation (HI), we obtain 



/(Ax) 



J P^'p(Ax;v' ,(J.')dfi' 



= < 



1 T^-^d/i' (Ax>0) 
A ■ 



A* 



e~^dfi' (Ax<0). 



(7) 



In this description, the renewal of // at the scattering and 
the translation to the following scattering point are treated 
as one step of the random walk. Although this p.d.f. have 
no information about fi', if we combine the information of 
the random walk with the scattering law, the information 
about (j, 1 can be taken out as described later. 

Corresponding to the p.d.f. f(Ax), the moment gener- 
ating function (m.g.f.) is defined as 



no) := 



c f(Ax)dAx 



Put 



1 + A* 



-djj! ', 



(8) 



which is the two-sided Laplace transform of the p.d.f. f(Ax) 
Clearly, /*(0) = 1. In general, at some values of 9 this in- 
tegral diverges. For example, for the large angle scattering 



model, the value of 9 



■jr- and 



gence points. However, in the following arguments, we need 
to consider only the value of 9 in the range, 

1 „ 1 



La 



< 9 < 



L- 



(9) 



and so we need not worry about such a divergence. 

We can derive the following formulas about the deriva- 
tives of f*(9) and the moments of f(Ax) from m.g.f. as 



d n r(9) 
d9 n 



(1-A*0) 



+ 



A 

P„'A* n 



(1 + A*< 



—dfi' \ , 



(10) 



Ax n := 



Ax n f{Ax)dAx = (-1) 



d9 n 



= n!|(-iry P^dn' + J P„,\ m dp! 

(11) 

In the range of 9 mentioned above (^|) , d d ^ > always holds 
and then /* (9) is a downward-convex function in this region. 
Therefore, the value of 9 — 9o 7^ at which f*(9o) = 1 
always exists. This quantity 9o, which is determined only by 
the functional form of f(Ax), includes the most important 
characteristics of the random walk in a moving medium. In 
fact, this 6>o together with v play an important role in the 
following description. Clearly, the sign of 9o is equal to the 
sign of v and Ax, i.e., it is positive throughout this paper. 



2.2.2 Random walk with absorbing barriers 

The problem of random walk in a moving medium with 
boundaries can be formulated as the problem of the random 
walk with absorbing barriers in the probability theory. The 
procedure is as follows. First, consider that a particle begins 
a random walk from the origin x = under the situation in 
which there are two absorbing barriers at the position x = a 
and x = —b (a, b > 0). Then the particle moves according to 
the p.d.f. f(x) at each step. If the particle advances beyond 
one of the barriers, it is absorbed by the barrier and the 
random walk is terminated then. 

To make the argument clear, we introduce the proba- 
bility density function of scattering points at the m-th step, 
f m (x). The integration of f m (x) outside of the boundaries 
(x < — b or x > a) denotes the probability that the parti- 
cle is absorbed at the m-th step. The integral of f m (x) over 
the inside region is not normalized to unity but denotes the 
probability that the particle is not absorbed until or at the 
m-th step and can carry out the m + 1 step. The function 
fm(x) can be written as a convolution of f m -i(x) by f(x) 
over the non- absorption region, 



f m (x) = / f{x — x')f m -\{x')dx. 



(12) 



Because the particle starts from x — 0, it can be seen easily 
fo(x) = 5(x) and fi{x) = f(x). All f m (x) can be calculated 
in principle using equation (113) iteratively. 

Next, we introduce the density of scattering points 
summed over all steps, 



are such diver- n{x) := ^ ' fm(x) {—b < x < a), 



(13) 



which constitutes one of the key concepts in the following 
description. If one is not discussing the number of steps be- 
fore absorption, the only relevant quantity is the density of 
scattering points n(x). We limit our considerations to those 
using n(x) in this paper, although the distribution of the 
number of steps before absorption is important to discuss 
the dispersion of acceleration time of particles. It is noted 
that n(x) is not the same as the physical number density 
of particles and their mutual relation is explained in the 
Appendix A. 

From equation ([l2]) and equation ([l3]), we can derive a 
Fredholm integral equation of the second kind for n(x): 
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n(x) = f(x) + / f(x — x')n(x')dz 

-b 



(14) 



Especially, for no-boundary case (a, b — > oo), the solution 
un(x) obeys 



tin(x) = f(x) + / f{x — x')riN{x')di 



(15) 



It is noted that n(x) satisfies certain integral conditions. 
To simplify the following description, we introduce the prob- 
ability density of absorption positions, 

h(X) = f(X) + J f(X - x)n(x)dx. (16) 

Here X instead of x is used for the spatial coordinate in 
the absorbing regions. n(X) satisfies the following two con- 
ditions. Since the particle is eventually absorbed by either 
of the absorbing barriers, 

-b poo 

h(X)dX + / ii(X)dX = 1. (17) 

30 J a 

In addition, in the problem of the random walk with absorb- 
ing barriers, a theorem which is called the Wald's identity 
holds. One of the simplest versions of this theorem is the 
following identity 

-b poo 

h(X)e- B ° x dX+ n(X)e~ e ° x dX = 1, (18) 

30 J a 

where Oo was defined in the previous subsection. In the shock 
acceleration, the barrier exists only at one side, i.e., a — > oo 
for downstream region and b — > oo for upstream region. 
In the following section, we use one of these conditions to 
normalize the approximate solutions. For a — > oo, we use 
equation (|l8|), i.e., 

n(X)e~ eoX dX = 1. (19) 

On the other hand, for b —* oo, we use equation (|l7j), i.e., 

n(X)dX = 1. (20) 



2.2.3 No-boundary solution and diffusion length 

In order to know the properties of random walk after large 
numbers of steps, it is helpful to consider the asymptotic 
form of un{x) for x — > +oo or x — > — oo. Using the con- 
volution theorem to equation (^), the two-sided Laplace 
transform of f m (x) for the no-boundary case turns out to 
be 

fW) = {ne)} m - (21) 

Therefore, from equation (|l3|), the two-sided Laplace trans- 
form of the scattering point density for no boundary case 
un(x) within the range < 8 < 8q becomes 



Thus, um{x) is calculated by the Bromwich integral as 



(22) 



n N (x) 



1 
27ri 



(J A- ZOO 



tin* (0)e 9x d9, 



(23) 



where a is real and < a < 6q. For x — > +oo, the integral 
path is taken on the —8 side and the residue at 8 — is 
dominant in the integral. The residue is calculated easily as 

a- 1 (8 = 0)= ' 



df{6) I 

de | e=o 



Then, we obtain the asymptotic form of njv (x) for x — » +oo 
as 



n N {x) 



df(B) 




de 





= =: no 

x 



for x — » +oo 



(24) 



The asymptotic form for x — » — oo is calculated in a simi- 
lar way. The integral path is taken on the +8 side and the 
residue at 8 = 8o is dominant in the integral. The residue is 
calculated as 



a-i(0 = 0o) = 



1 



a B Q X 



df'W I 

de | e=e 



Then, we obtain 

n N (x) = 



df*(6) I 
de \e=t 



&nx I Box 

-e u —: n e 



for x — > — oo . (25) 



For later convenience, we define the ratio of these two con- 
stants as 



no 
no 



df(e) 




de 






e=o 


df*(8) 




de 


e=e 



(26) 



If 8oL-,8oL+ 1 holds, p is almost unity (see equa- 
tion (10)). However, it is to be noted that p is not equal 
to unity, i.e., np ^ n' in general. 

Equation (123) indicates that the scale length of density 
of scattering points far upstream is equal to the inverse of 
do- We denotes this length as 



Ld - To- 



(27) 



While L+ and L_ are the scale lengths of single scattering, 
Ld is the scale length of multiple scatterings. If v is small 
enough and the scattering is isotropic, this scale length Ld 
should coincide with the diffusion length derived by the con- 
ventional diffusion equation, Ldo = k/V{ , where k is the dif- 
fusion coefficient (k — |An'), which is explicitly shown in 
Section |^. In contrast, when v is large, we can not safely 
use the conventional diffusion equation and the length I/do- 
However, even in such a case, we can regard Ld as a diffusion 
length, since the distribution is an outcome of multiple scat- 
terings in the random walk. Here and hereafter, we call this 
Ld simply 'the diffusion length' and call Ldo 'the conven- 
tional diffusion length'. We also show in Section ^ that this 
diffusion length Ld indeed coincides with the scale length of 
particle distribution derived by Peacock (1981) and Kirk & 
Schneider (1988) in the far upstream region for relativistic 
shocks. 

For later convenience, we also define the ratio of scale 
lengths of single scattering to the diffusion length as 



hi 

Id 



(28) 



a — toe 
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2.2.4 Probability density of pitch angle cosine at 
absorption 

As already mentioned, although the density of scattering 
points does not give direct information about pitch angle 
cosine of the particle fi' , we can derive the probability den- 
sity of fj! at absorption with an aid of the scattering law. The 
probability density of // at the absorption by the boundary 
x — a is given by 



e 1 * n(x)dx 



(29) 



P a (^) = ^e-^|l + J 

and that at the boundary x = — b is given by 

P b (ti') = P„>e~^ jl + J e-^n(x)dx^ 

The integral conditions for n(x) can be rewritten as the 
conditions for P a (p!) or Pt(fi'). For b — > 00, the condition 
( pp| ) can be rewritten in terms of P a (fjf) as 



(30) 



P a (tx')dfi = 1. 



(31) 



Similarly, for a — > 00, the Wald's identity ( |19| ) can be rewrit- 
ten as 



W) 



.A 



-dfj! — e 



-e b 



(32) 



Here, we define the total absorption probability by the 
boundary at —b 



Pab S (&) 



Pb(jj!)dp! 



(33) 



This function also means the probability that the particle, 
which is injected in downstream side at a distance of b from 
the boundary, reaches the boundary against the flow of scat- 
tering centres. If _R_ <C 1 holds, equation (B2J) also means 



Pabs(fe) = e 



(34) 



for any models. Since, for isotropic scattering and for v — » 0, 
Ld = 1 /0o approaches the conventional diffusion length Ldo 
as already mentioned, equation ( |34| ) agrees with the result of 
the conventional diffusion approximation obtained by Drury 
(1983). 



2.3 The return probability 

In order to apply the random walk theory described above 
to the shock acceleration, it is useful to consider the proba- 
bility density of the pitch angle cosine // when the particle 
which crossed the shock front with n' returns to the front 
crossing the reverse direction. For downstream, random walk 
with a finite b and a — » 00 is applied while for upstream that 
with a finite a and b — > 00 is applied. As will be shown in 
later sections, all the quantities associated with the shock 
acceleration can be represented with this probability den- 
sity. Hereafter, we denote the spatial coordinate by z with 
the shock front at z = in order to retain x coordinate 
for the random walk problem treated so far where x = 
corresponds to the injection point of the particle. 



First, let us consider the return probability for the up- 
stream. The particle first crosses the boundary toward (-) 
direction, i.e., from the downstream to the upstream with 
a pitch angle cosine fJ.'o(< —v). In order to calculate the 
probability density of the pitch angle cosine at the return 
jtt'(> —v), Ptjw(m'; Mo); it is useful to introduce the den- 
sity of scattering points before return, N-(z; /z ), which is 
defined on z < 0. This density can be calculated in the 
following way. The probability density of the first scattering 
point for the particle which crossed the boundary with pitch 
angle cosine fi is determined by the scattering law (^) . Par- 
ticles which have the first scattering point at z = —a make 
the scattering point density n(z + a; a, 00), where n(x;a,b) 
denotes the scattering point density of particles injected at 
x — in the random walk with absorbing barriers at x = a 
and x — — - b. Therefore, we obtain 

1 1 f°° 

N-(z; fi ) — — e A o + — / n(z + a; a, oo)e A o da (35) 
*a K Jo 

where z < and \q = X*(fj,' ). Using this, Puw(m';Mo) is 
calculated as 



Puw(At';Po) = Py. 



N-(z; fj!o)e^* dz 



(36) 



where A* = A*(^i'). For this case, the total return probability 
is always unity (equation ([jl|)), i.e., 

Puw(n'; n'o)dijf = 1. (37) 

We can rewrite this as an integral condition for N-(z; fJ,'o), 

1 r0 

P u , I N-(z; f^)e^ dzdn' = 1. (38) 



For the downstream, the probability density of pitch 
angle cosine fi'(< —v) at return for the particle which ini- 
tially crossed the shock front with fi'o(> —v) ((+) direction), 
Pdw(m'; Ak))i is calculated in a similar way. The scattering 
point density N+(z; n' ) is written as 



N+(z;n' ) = 4e A o + 1 



v o Jo 



n(z -b;oo,b)e x o db, (39) 



which is defined on z > 0. Using this density, Pr>w(fJ''', Mo) i s 
calculated as 



Pdw(m';Mo) = P^' / N+(z;n' )e ^ dz. 



(40) 



The total return probability -Pr(po) is n °t unity and is given 



^r(Mo) 



Pdw(m'; ^'o)d^' ■ 



(41) 



The integral condition (|32|) becomes a condition for 

N+(z;fi' ) as 

1 i-W N +l«M<-*W' = TTw' (42) 

as is shown in Appendix ^. 

In the following sections, we derive approximate ana- 
lytic solutions of Puw(m'i Mo) an d Pdw(//; (4i) where we use 
the normalization condition for the scattering point density 
iV_ and iV+ . 



© 0000 RAS, MNRAS 000, 000-000 



Random walk theory of shock acceleration 7 



3 APPROXIMATE SOLUTIONS 



where 



In this section, we derive approximate analytic solutions for 
the probability densities by obtaining approximate solutions 
for the integral equation (1l4f) . We start from the asymptotic 
solution of ( pi| ) and ( p^ ) for the no-boundary problem. First 
we approximate the no-boundary solution tin(x) as 



I n e ° x (x < 0) 



(43) 



Note that this approximation neglects a peak around x = 
arising from the injection (see Fig. |jj in Section ^| for Monte 
Carlo results). 

To obtain approximate solutions for problems with fi- 
nite boundaries from those for no-boundary problems, we 
transform the original integral equation ( p^| to an alterna- 
tive integral equation 

n(x) = un{x) — k(x, 0) — / k(x, x')n(x')dx' , (44) 
J -b 

where we define the kernel 

k(x,x') := / n N (x - X)f(X - x')dX 



+ 



n N {x- X)f(X -x')dX 



(45) 



(see Appendix |§| for its derivation). This integral equation 
is, of course, equivalent to equation (Q. 



3.1 Approximate solution of the probability 
densities 

First, we consider approximate solution for the upstream 
N- (z) taking b — > oo and a finite value of a. Starting from 
the approximation (^3|), the kernel of the integral equation 
(45) for n(x) becomes separable into x, x as 

k(x,x') =n' Q e Ba{x ~ a) (j ) -{x' -a) 

where 



<j>-(x'):= I f(X - x')e- e ° x dX. 



Substituting this into equation ([f4|), we obtain 

n(x) = n N (x) — n' o e 0o< - x ~ a ^ {(f)-(-a) + Ca} 

= n N (x) - n' o e 0o< - x ~ a) C' a 

where 

C a '■= I <t>- {x — a)n(x')dx' 



and 

C' a := <b-(-a) + C a . 

The term proportional to C' a represents the effects of the 
existence of the absorbing boundary at x = a. Although 
we do not need to explicitly evaluate C a in our procedures 
described below, C a can be solved if we substitute (46) into 
the definition of C a as 



C a 



1 



1 + n' C 2 
© 0000 RAS, MNRAS 000, 000-000 



<f>- (x — a)nN(x)dx — n C24>- (—a) 



C 2 := I <f>-.(x')e e ° x ' dx' . 



If we use the approximation (^3J) again and remember- 
ing Ld = l/#o, the equation d35|) becomes 



N_(z;v' ) = C (AS)^e A o +C 1 (AS)-^e 1 



where 



C (\*o) := 1 + n \* - 



(46) 



(47) 



and 



Ci(Aj) := 



n' /6o «o_ 

1 - \q 6 \* 



C' a {a) 



x o da. 



In equation (^), the second term with the scale length Ld 
represents the distribution formed by the diffusion process 
after a large number of scattering, similar to the asymptotic 
form of un(x) for far upstream (|25[). The first term repre- 
sents the correction to the diffusive term by the following 
three non-diffusive effects; the density of scattering points 
contributed by particles which take only a few steps, the ef- 
fect of finite initial mean free path Aq (finite distance from 
the boundary to the first scattering point), and the absorb- 
ing effect near the boundary. When we take the limit AS — » 0, 
this term disappears and equation (^) is just reduced to the 
diffusive solution. Depending on AS, Co can take a positive 
or negative value. For small Aq the correction term is pos- 
itive and it mainly reflect the scattering points of the few 
step particles, while for large Aq it is negative accounting 
for the effects of finite mean free path and absorption at the 
boundary. These features are mentioned again in Section ^ 
for a specific scattering model. 

Substituting this to equation (|3(j) , we obtain an approx- 
imate solution for Puw 

Pvw^'^'o) = Co(Ag)^l + Ci(Ag) ,f M ' A * . (48) 



A* + A ( * 



A* + Ld 



However, this approximate solution does not satisfy the con- 
dition ( j37|) in general because we do not use the exact so- 
lution of n(x). Although there may be various ways of rem- 
edying this problem, here we adopt a simple and practical 
way to renormalize the amplitude of Ci so as to satisfy j38| ) 
for a given Aj$. Thus, instead of Ci we use the coefficient Ci 
defined by 



Ci(AS) = 



1 - C 05+ (AS) 



g+(L-D 

where we define 



dfj,' . 



(49) 



(50) 



This method also has a merit in avoiding to calculate a com- 
plicated integral in the original definition of C\. 

Next, we consider approximate solution for the down- 
stream (z) in a similar way by taking a — > oo and a finite 
value of b. Starting from the approximation (ji^), the kernel 
becomes 

k(x,x') = no4>+(x') 
where 
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<j>+{x) := / f(X-x)dX. 

J — oo 

n(x) is obtained as 
n(x) = tin(x) — noC' b 
where 



(51) 



C' b := <j>+(0) + I <j>+(x')n(x')dx' 



Substituting this into equation (^9|) and using the approxi- 
mation (^) again, we obtain 



N+{z;n' )=C' ^e ^ + C[ 
A o 

where 



C'o(K) := 1 - n Xo + ■ 



1 + 6oK 

C[(\* ) :=n { 1- / C£(6)e A S rf & 



(52) 



(53) 



Similarly to the case of N- , the second term in (£32|) denotes 
the diffusive term, which corresponds to the asymptotic form 
for far downstream (^4|), and the first term represents the 
correction to it. C' takes positive value for small A ( * and 
negative value for large Aj> Using this to equation (MCJ) , we 
obtain approximate solution of Pdw 

Pbw(//; fi' ) = C (X* ) y + Cfl(Ag)iVA'. (54) 

A + A ( , 

We again renormalize C[ so as to satisfy the integral con- 
dition ( J42] ) for the same reason as in the upstream. Thus, 
instead of C[ we use C[ defined by 



g+{L D )(X* + Lr>) 
where we define 

Pu'X* , , 

-dfx 



(55) 



(56) 



and we use the relation 

- ff _(-L D ) = ff+ (L D ), (57) 

which is derived by equation (8) together with the definition 
of 6o and J P^tdfx = 1. The total return probability Pr 
is calculated by equation ( pi] ) easily, 

Pr(mo) = CS(AS)p-(AS) + Ci(AS)h- 
where we define 



(58) 
(59) 



h- := J P^X*dfjf. 
3.2 Multi-step approximation 

To compare our approximation with more conventional dif- 
fusive approximation, for reference we consider also the lat- 
ter approximation which we call multi-step approximation. 
This is done by ignoring the first term in equation (^) or 
([52]), i.e., by setting Co, Co — > 0. In this approximation, the 
density profile takes the same form to the results of the dif- 
fusive approximation. But, since we do not take the limit 



A2 — > 0, this approximation includes partly the effect of the 
finite mean free path. (Although, the effect of the first step 
return is lost.) 

Thus, in the multi-step approximation, the probability 
density functions become 

Puw(m';Mo) = — 77 — 7TT7~T — ' ( 60 ) 



g+(Lu) A* + L D ' 



Km'; 



Pu'X* 



Mo) 



g+(L u ) Xq" + L u 



and 



PR(tA) = A 



Ag + Ld ' 



where we define 
A 



(61) 



(62) 



(63) 



In the multi-step approximation, the pitch angle distribution 
at return becomes independent of no and n . Further Puw is 
independent of fi' an d Pdw is separable with respect to fi' Q 
and p! , Although in general the multi-step approximation 
is not satisfactory, there are cases where this makes some 
sense as shown in Section m 



3.3 The absorption probability 

The absorption probabilities of the particle injected at a dis- 
tance b in the downstream from the shock front, Pb{(J>') and 
Pabs(f>), are of some theoretical interest, although it is not 
explicitly used in later considerations. To obtain their ap- 
proximate expression, we start with equation ( p"l| ) together 
with the approximation (^) . We renormalize the amplitude 
of c := —noC' b so as to satisfy the Wald's identity (B2) as 



g+(L n ) I L D 



1 --^ 

e d 



JVC (A')e-s*_ d , 



Ld — A* 



* [-" P„,A» , 



(64) 



Using c, we obtain 



at//) = />„ ^ e-&C (X*)+n' e 4 + £A* }. .iOM 



P abs (6) = 



P M ,C (A*)e- — <V 



+Ld3+(£d) |n e ^ + cAj 



(66) 



It is to be noted that P a bs(6) has the scale length Ld for 
large 6. 

The absorption probability of the particle injected in 
the upstream at the distance of a from the shock front 
and the pitch angle cosine fjf , P a (//), can be also calcu- 
lated in a similar manner. Using equation (46) with ( pf3[ ) 
we renormalize c' := n C' (a)e~ 6 ° a so as to satisfy equa- 
tion (H), although we do not present the result here. In 
principle it is possible to calculate Puw and Pdw by using 
these P b (fi'), P a (/x'). However, since the expressions are too 
cumbersome and impractical to denote Puw and Pdw in an 
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analytical form, we use the equation fl48p and (|54|) for Puw 
and Pdw respectively in the following argument. 

We make further approximation such as the multi-step 
approximation. If we set un{x) = no for all x, which is the 
diffusive solution, and if we ignore the first term in equa- 
tion @, using Wald's identity (^), we obtain 

P b ^')=Ae-^^-, (67) 



Pabs(fe) = Ae 



(68) 



In the limit of v — > 0, in which the conventional diffusion 
approximation holds, the constant A approaches unity and 
this expression of P a bs(b) agrees with equation (Q). Pea- 
cock(1981) derived the absorption probability in a similar 
form (see his equation (24); he did not give a concrete ex- 
pression of the constant A), but did it in the downstream 
fluid frame under non-relativistic fluid speed condition and 
the diffusion approximation. In the following section, we will 
make comparisons with his expression for the isotropic scat- 
tering. 

Corresponding approximation in the upstream can be 
made in a similar way. If we set the diffusive solution 
un{x) oc e L o , n(x) becomes c'e L o , where c' is some con- 
stant. Ignoring the first term in equation (^) and determin- 
ing c' through equation (pll), we obtain 



P„j\* 



g+(L D ) A* + Ld 



(69) 



If we use these expressions for Pt(//) and P a {^') to calculate 
Pdw and Puw, the results coincide with the results of the 
multi-step approximation ( pu) ) , (^lj) obtained in the previous 
subsection as should be since only diffusive return is taken 
into account. Therefore, we call these expressions for Pf, (//'), 
Pabs(&) and P a (p.') 'the multi-step approximation', too. 



4 RESULTS FOR ISOTROPIC LARGE ANGLE 
SCATTERING MODEL 

In this section, we apply the random walk theory developed 
so far to a specific model of scattering. We consider here the 
model in which the mean free time r is independent of /i' 
(although the energy dependence may be allowed, it is not 
relevant here), i.e., t(/j,',v') — To(v'), and the scattering is 
isotropic in the fluid frame, 



A — v \fi jr , 



P / = - 
2' 



(70) 



This model has been widely used both in analytic works 
(Peacock 1981) and in the Monte Carlo simulations (Ellison, 
Jones & Reynolds 1990) . To check the validity of our model, 
we also perform the Monte Carlo simulation in which the 
particle position is traced step by step faithfully according to 
the adopted scattering law. In the simulation, the escaping 
boundary is set at X csc — 15Ld, which is far enough not 
to influence the return probability of the particle (see e.g. 
equation (p8|)). 



4.1 Properties of the random walk with isotropic 
scattering 



In this model, A* is given by 

a* = r f t/To|// + u\. 



(71) 



For the reason mentioned in section 2.1, we can take the 
unit of length as Tfv'ro in the following argument. Thus, A* 
becomes 



A* = + v\ 
and 

L+ = l + f, L_ =1 — 1/. 

The p.d.f. of the random walk is given by 

l El (^) (*>0) 

= < 1 \x\ 

±Ex(f-) (x<0) 



(72) 



(73) 



(74) 



where E n (x) is the ex pone ntial integral defined by (see 
Abramowitz & Stegun (1965)) 



Ejx) 



t*> 



-dt. 



The m.g.f. is given by 

f {d) -26 ln {i-e(i-v) 



(75) 



(76) 



1 _i_ 2 

3+^' 



—2 

■ X 



(77) 



and the moments of f(x) are calculated as 

1 

3' 

Fig. shows the p.d.f. f(x) (a), and the m.g.f. f(6) (b) 
for v = 0.1,0.5,0.9. It is seen that the p.d.f. becomes more 
asymmetric for larger v because of the advection effect. 

The equation to determine #o (or Lo := l/#o) is given 

by 



26»o - In 



l+fl (l+i 

1 - 0o(l - i 



0. 



(78) 



Fig. ^ presents the scale lengths of the random walk, L+, 
L_ and Ld- We also show the conventional diffusion length 
Ldo mentioned in Section 2.2.3 extrapolating it to high v. 
Since it depends on the Lorentz factor of fluid speed Tf in 
the unit of length used here (Ldo = l/(3Tiu)), we show the 
two cases; the case of a relativistic particle in a relativistic 
flow (v' ~ 1, v ~ Vf \ lower dotted curve), and the case of a 
non-relativistic particle in a non-relativistic flow (v',Vf -C 
l,Tf ~ 1; upper dotted curve, Ldo = 1/(3^)). The diffusive 
scale length for a pitch angle diffusion model Lpd (Kirk & 
Schneider 1987a) is represented, too. (Here, we choose the 
pitch angle diffusion coefficient in their model so that it has 
the same spatial diffsion coefficient as our model.) For v 1, 
Ld can be approximated as 

L D ~ — , (6»o ~ 3u) (79) 



3i/ ' 

and for v — * 1, 



L_(l 



1 — v ' 



(80) 



where 
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Figure 1. (a) The p.d.f. f(x) and (b) the m.g.f. f*(0) for isotropic large angle scattering model. Those for v = 0.1 (solid curve), 0.5 
(dashed curve) and 0.9 (dotted curve) are shown. 



c 

CD 




Figure 2. The scale lengths of the random walk for the isotropic 
large angle scattering model in the unit Y^v'tq. The diffusion 
length Lo = l/0o (solid curve), the scattering length to the up- 
stream direction L- (lower dashed curve) and that for down- 
stream direction L^. (upper dashed curve) are shown as functions 
of v, the ratio of the fluid velocity to the particle velocity. The con- 
ventional diffusion length l/(3Ffi^) (y = Vf, lower dotted curve) 
and 1/(31/) (r"f = 1, upper dotted curve) are also shown (see text). 
Lpo (dot-dashed curve) is the diffusive scale length for a model 
of pitch angle diffusion derived by Kirk & Schneider (1987a). 



e: =ZT e 



(81) 



Therefore, if the ratio v is small enough (y 0.3), our dif- 
fusion length Lo coincides with the conventional one Loo, 



as expected. On the other hand, if v becomes larger, the 
difference between Lo and Loo becomes larger. As men- 
tioned in Section 2.2.3, for large v, Loo becomes meaning- 
less because the conventional diffusion approximation is no 
longer valid. However, Lo still has a meaning of 'the diffu- 
sion length' and this is explicitly shown later in Fig. Rln 
fact, this scale length has been derived by Peacock (1981) 
and Kirk & Schneider (1988) independently as the scale 
length of the particle distribution on the far upstream re- 
gion of relativistic p arallel shocks for large angle scattering 
models. In Peacock (1981), the equation that is equivalent 
to equation ( f78| ) was derived from the conservation law of 
particles in the upstream of the shock using 'the relativis- 
tic diffusio n appr oximation' (which was named by Kirk & 
Schneider (1987a)), where (a + m)/(l + u\) corresponds to 
our 9o in his equation (16). In Kirk & Schneider (1988), such 
equivalent equation to Q7q) was also derived as the equation 
to determine the eigenvalues for large-angle scattering op- 
erator. In their equation (34), corresponds to 8o and lo 
was taken to unity. The upstream density profile obtained 
in a previous Monte Carlo work for oblique shocks (Naito & 
Takahara 1995) deviates from the diffusion approximation, 
which can be interpreted using Lo obtained here if we con- 
sider in the de Hoffmann- Teller frame. Fig. ^| also shows that, 
when v becomes larger, the difference between L+ and L_ is 
larger and Lo approaches to L_. Behaviour of Lpd is simi- 
lar to that of Lo but takes a smaller value. It indicates that 
the diffusive scale length depends on the scattering model 
explicitly. 

The constants no and n' Q , which determine the density 
of scattering points in the no-boundary case, becomes 



1 



no 



"o 



1 



- 2v 



■0o. 



(82) 



6»o(l-^ 2 , 

Fig. [] shows the ratio p — n' /no as a function of v. For the 
limit v — » 0, p approaches to unity, which agrees with the 
results of the conventional diffusion approximation (Drury 
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Figure 3. The ratio n' Q to no as a function of v. 




Figure 4. The density of the scattering points for no-boundary 
solution rijf{x) (solid histogram) for v = 0.5 obtained by 
a Monte Carlo simulation and the asymptotic approximation 
n^p (x) (dashed curve). 



1983). However, as v becomes larger, p becomes smaller and 
approaches as 2e/L 2 _ when v — > 1. 

Fig. ^ shows an example of the no-boundary solution 
un(x) (solid histogram) for v = 0.5, which is obtained 
by the Monte Carlo simulation, together with the asymp- 
totic approximation n$(x) (equation (pL3|); dashed curve). 
It confirms that n^r (x) approaches njv (x) asymptotically for 
x — » — oo or x — » +oo. Although tin(x) has a peak around 
x — 0, our approximation ignores this peak as described in 
Section ^. Fig. |H] shows an example of the results when the 
single boundary is located at a = 5Lo(b — > oo). The density 



Figure 5. The density of the scattering points n(x) (solid his- 
togram) obtained by a Monte Carlo simulation is shown for 
v = Vf = 0.5, a = 5Lj),b — > oo together with our approximate 
solution n(°) (x) and that of the conventional diffusion equation 
n^iff(x) (dotted curve). 



of scattering points n(x) (solid histogram) for v = 0.5 is 
obtained by a Monte Carlo simulation, along with our ap- 
proximate solutions and that of the conventional diffusion 
equation (Drury 1983; his equation (3.41)) ndiff(a;) (dotted 
curve). It is clearly seen that the conventional diffusion equa- 
tion does not give a sufficiently good approximation to the 
real solution for high v. Especially, the random walk result 
shows that n(x) has a finite value at the absorption barriers, 
while ndiff (a;) becomes there. The small difference between 
the Monte Carlo result and the our approximation near the 
absorbing boundary is mainly caused by omitting the peak 
in njv(x) around x — in making our approximation. 



4.2 Return probability densities 

Here, we present the return probabilities of particles for the 
isotropic large angle scattering model under the approxima- 
tion made in the previous section. This model gives 



h- : 
and 

A = 



4 



l 



where we define 

, , . x — ln(l + x) 
4>(x) :=2 \ i. 



(83) 



(84) 



(85) 



(86) 



(87) 
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Figure 7. The absorption probability from z, P a bs( 2 )i f° r v = 
0.7. Results of our approximation (solid curve), the multi-step ap- 
proximation (dashed curve), and Monte Carlo simulation (filled 
circle) are represented. The results using the length l/(3Tfu) 
(lower dotted curve) or 1/(3^) (upper dotted curve) instead of 
Ld in the multi-step approximation are also shown as in Fig. 2. 
The curve 1 / (3u) also corresponds to the extrapolation of Pea- 
cock's expression (his cq. (24)) to relativistic fluid speed. 



From these quantities, the probability densities of pitch an- 
gle at return Ptjw and Pdw are calculated by equation (Jii]) 
and j&4]), while those for the multi-step approximation are 
calculated by equation ( |6c| ) and (JsTJ) , respectively. The total 
return probability Pr is also calculated by equation or 
equation (|62|) . 

Fig. Hshows the probability density of pitch angle at 
return Ptjw(mo) m ( a ) an d Paw(nb) in (b) for v = 0.5 
and for several values of the initial pitch angle cosine fj,' . 
Our approximation (solid curves) gives a fairly good fit to 
the Monte Carlo simulation (grey histograms) for all initial 
pitch angle cosine fi'o. For small /xo, the fit deviates a lit- 
tle because the effects of the peak in un{x) around x = 0, 
which is ignored in our approximation, affect the probabil- 
ity density. The multi-step approximation (dotted curves) 
does not fit well when the initial mean free path Aq is small 
because effects of the return after only a few steps of scat- 
tering becomes important. But, for large Ao and for Pdw, 
the multi-step approximation can fit to the Monte Carlo 
simulation. 

The expression of the absorption probability P a b s (&) by 
our approximation is a little cumbersome and is presented 
in Appendix The absorption probability by the multi- 
step approximation is given by equation ( |68| ) together with 
equation (pif). Fig. ^ presents P a bs(&) as a function of the 
initial position 6 for v = 0.7. The result of our approxima- 
tion (solid curve) gives a good fit to the Monte Carlo results 
(filled circles) (small difference is in part caused by omitting 
the peak in tin(x)). The multi-step approximation (dashed 
curve) deviates from the Monte Carlo results, but gives a 
correct scale length (L-d). This result directly confirms that 



Figure 8. The total return probability Pr(/^o) for v = 1/3 from 
the downstream. Here, the initial pitch angle cosine fio is mea- 
sured in the boundary rest frame. Results of our approximation 
(solid curve), the multi-step approximation (dashed curve), Pea- 
cock's approximation (dotted curve) and Monte Carlo simulation 
(filled circle) are shown. 



the diffusion length for such large v is given by our diffusion 
length Ld. Results using the conventional diffusion length 
Loo instead of Ld in the expression of the multi-step approx- 
imation (equation (^) and (|6^)) for two case as in Fig. g are 
also represented; lower dotted curve for l/(3T{is) and upper 
dotted curve for 1/(3^)). It should be noted that if we ex- 
trapolate the expression (24) in Peacock (1981), which was 
derived by adopting the diffusion approximation for parti- 
cles in downstream of the shock front, to relativistic fluid 
speed in fluid frame and then transform it to the expres- 
sion for the boundary rest frame, it turns out to correspond 
to the multi-step approximation ( p8| ) except that it has the 
scale length of return l/(3v) (in our unit). These results do 
not give a satisfactory fit. 

Fig. ^presents the total return probability Pk(^o) from 
the downstream for mildly high v (y — 1/3) as a func- 
tion of initial pitch angle cosine no, which is measured in 
the boundary rest frame. The results of our approximation 
(solid curve) excellently agree with the Monte Carlo simula- 
tion (filled circles). The do tted curve shows the approximate 
solution of Peacock (1981) (shown in his fig. 2). The result of 



the multi-step approximation (dashed curve) , which uses the 
correct length Ld, is similar to Peacock's one, but slightly 
better than it. These two diffusive approximations do not 
agree with the Monte Carlo results for small ^o- The reason 
is similar to that for Ptjw and Pdw; in diffusive approxima- 
tions, effects of return after only a few steps of scattering 
are not accounted for. In general, for higher v or smaller Ao, 
such effects become important and these diffusive approxi- 
mations become worse. 
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Returning Pitch Angle \i' Returning Pitch Angle n' 

Figure 6. The probability densities of pitch angle at return for u = 0.5 (a) Puw (m' ; Mo ) an d (b) Pdw (a 4 ' i Mo ) f° r severa l values of 
the initial pitch angle ^t Q . Results of our approximation (solid curves), the multi-step approximation (dotted curves) and Monte Carlo 
simulation (grey histograms) are shown. 




4.3 Density of scattering points 

Here, we show the results of the density of the scattering 
points for the particles which cross the shock front with the 
pitch angle cosine fi' . Fig. ^ shows N-(z;fi' ) in (a) and 
N+(z; Hq) in (b) calculated by the Monte Carlo simulation 
(histograms) for several values of the initial pitch angle n' 
along with our approximate solutions. It is seen that our ap- 
proximation (solid curves) coincides well with Monte Carlo 
simulation. But, a discrepancy is recognized for z > 0.3L_ 
in (b). Because the normalization conditions used in our 



approximation (the equation (p7|) , (J42|) ) employ only the 
information of finally absorbed particles and do not care 
about escaping particles to the far downstream, our approx- 
imation fails to give a good fit in a distant region where 
escaping particles dominate, while it gives a satisfactory fit 
in a near region where almost all the last scattering points 
exist. This also means that our approximation works well for 
absorbed particles, therefore the return probability densities 
themselves are represented by our approximations quite well 
as shown in Fig. n and Fig. H. 

Another important feature of N+ and iV_ is that the 
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density of the scattering points near the boundary depends 
on the initial pitch angle fi' . Here, we explain this feature 
for N+ (that for iV_ is similar). Fig. | (b) shows that N+(z) 
increases when z — > for small A2 but that it decreases when 
z — > for large A ( *. This can be explained as already men- 
tioned in Section 3; for small Aq the scattering point of few 
step particles are important, while for large Aq the effect of 
finite initial mean free path is important. In our approximate 
expression of (p^), while C[ is always positive, C' becomes 
negative for large Aq, as shown in Fig. ho|(b) . (This feature 
also appears for N-, as shown in Fig. fL0| (a).) Because of 
this feature, when we average the return probability over an 
incoming particle distribution, the first term in (|52|), which 
denotes the non-diffusive effects, will become small enough 
compared with the diffusive term which is associated with 
C[. Then, there appears a possibility to use the multi-step 
approximation to determine the averaged pitch angle distri- 
bution over all incident particles even if u is large. The same 
consideration applies for N-, too. As will be shown in the 
following section, at least for the isotropic large angle scat- 
tering model and for typical incoming fi distributions, this 
is actually the case. In practice, for this reason, we can use 
the multi-step approximation to determine the pitch angle 
distribution at the shock front for arbitrary shock speed as 
shown in the following section. 



5 APPLICATION TO THE SHOCK 
ACCELERATION 

In this section, we apply the random walk theory developed 
so far to the shock acceleration and calculate the spectral 
index of the accelerated particles for the isotropic large an- 
gle scattering model described in the previous section. In 
the following, we consider the acceleration in parallel shocks 
where the fluid speed measured in the shock rest frame are 
uniform on each side of the shock front; V u for the upstream 



and Vd for downstream respectively. We also consider only 
relativistic particles («' ~ 1, i.e., v ~ Vf) because, in this 
case, the ratio v does not change on each side of the shock 
front even if the particle gains energy in the process of multi- 
ple crossing at the shock front. In addition, for such particles, 
the Lorentz transformation of /i to /i' becomes as simple as 



M- Vt 
1-ltVt 



dfj,' 



Vi 1 



(i - my 



dp, 



(88) 



where fi is the pitch angle cosine measured in the shock rest 
frame. Hereafter, we use subscripts u, d or superscripts (u), 
(d) to indicate quantities in the upstream and downstream 
regions, respectively. The pitch angle cosine measured in the 
upstream fluid frame is denoted by ^S u ' and that measured 
in the downstream fluid frame by fi^ . 

In the following, we present the results only for two 
simple jump conditions r = V u /Va = 4 and r = 3V^ 2 
(i.e., V u Vd = 1/3). But, since the results are expressed in 
terms of V u and Vd, they are applicable to any jump condi- 
tions. (Jump conditions for re lativistic s hocks are discussed 
in Blandf ord fc McKee fll976| ), Peacock ([t98l|) , Appl & Ca- 
menzind (|l988|) and Heavens & Drury ([19881).) 



5.1 Method of determining the pitch angle 
distribution at the shock front 

The spectral index of the accelerated particles a is given by 



(In 5) 



(89) 



where (Pr) is the averaged return probability for the parti- 
cles crossing the shock front toward downstream and (In 5} is 
the logarithm of the energy gain factor (see Peacock 1981). 
(Here, () denotes to average over the pitch angle distribution 
crossing the shock front.) (In 8) is given as 
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(In 6) = 2 In F rcl + (In 5) ud + (In 8) du , (90) 
where 

(lnS) ud := ^(l + KoiM^)), (91) 

(ln«5) du := (ln(l - KdMdu)) (92) 

and Kei is the relative velocity of the upstream fluid with 
respect to the downstream fluid, 

v u -v d 



v iel = 



VuV d ' 



r r( ,i :— 



\/l - Kc 



(93) 



If the dependences on [i and v' in A and P^i are separable, 
these quantities are independent of the particle energy, i.e., 
v . 

Introducing the normalized pitch angle distribution of 
particles crossing the shock front from upstream toward 
downstream by $ u d (//"'), and that from downstream to- 
ward upstream by $du(A £< ' d ' ) )> i- e -> 



$ud(M (u) )V u) = / $ du (M (d) )V d) = 1, 



we can write relevant quantities as follows. Using (p^) 
and <I > ud (/i' d - ) ), (Pr) is written as 



(p R ) = / p R d V d) )* ud (M (d W d) . 



(94) 



The logarithm of the energy gain factor is given by 



(ln<5) ud = I ln(l + V Icl ^) "H /; ^ud; #ud(/x (u) )dM (u) (95) 



and 
(In S)du 



(96) 



The factor Pj^\^ d ) / (Pr) in equation ( |95| ) means that only 
particles which return again to the shock front are to be 
considered. 

The pitch angle distribution $ud(M <u> ) an d <l? d u(/z^) 
in a steady state are determined as follows. If $ u d(/i ) is 
known, $du(£t' d ^) is determined by & u d(fJ>^), 



"I>du(M 



1 p(d) / (d). ,,(dV 



4 d) (/4 d) ) 



(97) 



On the other hand, Qud(l-i^) is also determined by 

*du(M (d) ), 



$ud(M (U) ) = 



P^0^o)$ du (/^)<K 



(d) w „(d) 



(98) 
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Figure 11. The kernel of integral equation (100) for V B = 0.7 
and r = 3V^. The numerical evaluation based on our approxima- 
tion for fiQ = —0.01 (solid curve), —0.5 (dashed curve) and —1.0 
(dotted curve) are shown together with the result of the multi- 
step approximation (thin dotted curve). The results turn out to 
be almost independent of fio and all curves give similar results. 



This kernel denotes the probability density of the particle 
which crosses the shock front from downstream to upstream 
with /Xq crosses the shock front again from downstream 
to upstream with /i' d ' after one cycle. If $du(/-t ) can be 
solved, $ud(A t ' u ') is calculated through equation (pi]). Al- 
though in the above we have made the integral equation for 
$du(£i' d ^), one can alternatively make the integral equation 
for $ ud (A* (u) ). 



5.2 Approximate solution of the pitch angle 
distribution at the shock front 

It is generally difficult to solve the integral equation ( |99[ ) 
analytically even if our approximation are used. Then, we 
numerically evaluate the kernel for our approximation. On 
the other hand, in the multi-step approximation, <& ud ([i) and 
$du(M) can be calculated analytically without solving the in- 
tegral equation (|99|). Because the dependence of PdwKm'j Mo) 
in equation ( pil] ) are separable in fjf and fi' in this approxi- 
mation, it becomes 



$ du (M 



(dh 



P.'Xd 



(101) 



Therefore, in a steady state, $du(A t ^ d ' ) ) satisfies an integral 
equation (a Fredholm equation of second kind) 

* d u(M (d) ) = J " d ^(M (d) ;^ d) )<i>du^ d V^ d) , (99) 

where the kernel of this integral equation is given by 
*(M Cd U d) ):= f { ( C'£ ^ OMi U U u y/4 u) .(100) 



Because Puw(m'; Mo) in equation ( |6o| ) is independent of fjf 0) 
we obtain 



<Mm (u) ) 



5+ (L D )A*+L D - (1 ° 2) 

These results agree with the results obtained by Pea- 
cock (1981}, but Ld was replaced by Ldo in equation (101). 

Fig. [ll| shows the kernel of equation (^) for V s = 0.7 
and r — 3V S 2 . The numerical evaluation based on our ap- 
proximation for no = —0.01 (solid curve), —0.5 (dashed 
curve) and —1.0 (dotted curve) together with that of the 
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multi-step approximation are shown. It is seen that the ker- 
nel is almost independent of the initial pitch angle cosine 
Ho- The reason for this coincidence is considered to be due 
to compensation effects of several non-diffusive effects when 
averaged over the initial pitch angle as was mentioned in 
the previous section. Although this is partly due to our as- 
sumption of the isotropic large angle scattering, it is of a 
great importance in applications to the shock acceleration. 
Although the multi-step approximation is a poor approxi- 
mation to the detailed particle transport, it can be a good 
approximation when we average over the pitch angle distri- 
bution. Thus, this result gives a justification for using the 
multi-step approximation even for relativistic shocks. 

Fig. [12J presents the pitch angle distribution at the shock 
front in various frames for various shock velocities. The re- 
sults of the multi-step approximation (dotted curve) and the 
Monte Carlo simulation (solid histogram) are shown. As is 
seen, the results of the multi-step approximation agree quite 
well with the Monte Carlo results even for highly relativistic 
shocks. This is not surprising once we admit that the ker- 
nel is well approximated by the multi-step approximation. 
Thus, we have confirmed the point mentioned above that 
the multi-step approximation gives a good fit to the steady 
state pitch angle distribution at the shock for any shock 
speed. Also, it is to be noted that the assumptions used in 
Peacock (1981) are adequate for arbitrary shock speed if the 
correct diffusion length Ld is used in the downstream as in 
the multi-step approximation. 




2 3 4 5 6 7 

Compression Ratio r 



Figure 14. The spectral index as functions of the compression 
ratio r for various shock speed V B . The solid curves are our results 
based on the multi-step approximation, the dotted curves are the 
results of fitting formula of Ellison et al. (1990) and the filled 
circles are our Monte Carlo results. The dashed curve is the results 
of the conventional diffusion approximation, i.e., a = (r + 2)/(r — 
I)- 



5.3 Approximate solution of the spectral index 

In the previous subsection, we have shown that the pitch 
angle distributions at the shock crossing $ u< j(At) and $du(M) 
can be well approximated by the multi-step approximation. 
Here, using this approximation, we show analytical expres- 
sions of the quantities associated with the shock accelera- 
tion, {Pr), (In 5) and a, for the isotropic large angle scatter- 
ing. First, using Prt(/io; Vd) (equation (|u2[)) and $ u( j(/i) by 
the multi-step approximation, (Pr) is calculated as 



(P R ) = A w -U(V u ,V d ), 



where 
U(V u ,V d ) := 1 + 

with 

Xu := 



Xd 



Xu ~ Xd 



r u 2 (v u -v d )L[ 



Xd := f + 



r d 2 (v u -v d )4 d) ' 



Xo := 



Vu-Vd 



Xi := 



Vu-Vd 
l + V d ' 



(103) 
(104) 

(105) 

(106) 

(107) 
(108) 



R ■= xoXd (R+ = xoXu)- 

i/j is defined in (|S7|), 

Then, the logarithm of the energy gain factor is given 

by 



fe-5^)VK%o) + (Xd-l)*(Xd)-Uu-l)*(Xu) 
-L Aid Xu/ (109) 

(xd-i)*(«)-(xu-i)V'(-RV ) 

and 

(ln5>du = ln(l + Kei) + ^(Xi) - \, (HO) 
where 

ln(l + xo) In(z - 1) + dilog (^) - dilog ) 

with Dilogarithm dilog (x) defined by (see Abramowitz & 
Stegun 1965) 



dilog (a;) 



lnt 
t - 1 



dt. 



(112) 



(In 5) ud 



ln(l + Kc 



Using these results, ( p9| ) gives the spectral index a. In 
Fig. |l| (a) are presented the quantities associated with the 
shock accelerations as functions of shock speed V s for the 
compression ratio r = 4. In Fig. [l3| (b) are shown the re- 
sults for different jump condition (KiVd = 1/3). As is seen 
in these two figures, the result of the multi-step approxima- 
tion excellently agree with the results of the Monte Carlo 
simulation even if shock speed becomes highly relativistic. 

Fig. [li] presents the spectral index a as functions of the 
compression ratio r similar to fig. 4 of Ellison et al. (1990). 
Our results (solid curves) basically coincide with their ap- 
proximate expression (dashed curves; their equation(29)). 
Small discrepancy between these two results may be caused 
by the difference of the method to determine the spectral 
index a; we calculate a using equation (^9|), while a was de- 
termined by fitti ng the calculated particle spectrum directly 
in Ellison et al. jl990| ). 
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[3.) in the Fluid Frame in the Shock Frame (Dj in the Fluid Frame in the Shock Frame 




Figure 12. The pitch angle distribution of the particles crossing the shock front. Figure (a) refers to the particles crossing from upstream 
to downstream while figure (b) does those from downstream to upstream. The pitch angle distribution both in the fluid frame and the 
shock rest frame are shown. The results of the multi-step approximation (dotted curve) and the Monte Carlo simulation (solid histogram) 
are shown. 




Figure 13. The quantities associated with the shock acceleration, the spectral index cr, return probability (Pr) from the downstream 
and the logarithm of energy gain factor (In 8) for the compression ratio r = 4 (a) and for the relativistic jump condition, V u Vd = 1/3 
(b). Curves show the results of the multi-step approximation, while filled symbols show those of Monte Carlo simulation. 
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6 DISCUSSION 

Although we have considered the particle acceleration in 
parallel shocks with the large angle scattering model, it is 
of some interest to discuss whether our formulation can be 
extended to other scattering regime. 

Recently, the acceleration in ultra-relativistic shocks 
has attracted some attention in relation to the ultra-high- 
energy cosmic rays (UHECR). However, in such highly rel- 
ativistic shocks, the large angle scattering model adopted 
here can not be used safely in a physical sense; because the 
residence time of particles in upstream is quite short (it is 
shorter than the gyro-period of particles) , its point-like fea- 
ture may lead to unsuitable results (Bednarz & Ostrowski 
1996). This is also problematic for the pitch angle diffusion 
model based on the quasi-linear theory; because it assumes 
the resonance between the magnetohydrodynamical waves 
and particle gyro-motion, it requires time much longer than 
the gyro-period. Gallant & Achterberg (1999) proposed two 
deflection mechanisms of particles in upstream; regular de- 
flection owing to the gyro-motion in uniform magnetic field, 
and 'direction angle' diffusion by randomly oriented small 
magnetic cells, where the direction angle is defined as the 
angle between the particle velocity and the shock normal. 
Gallant, Achterberg & Kirk (1999) and Kirk et al. (2000) 
calculated the spectral index using such direction angle dif- 
fusion model. Bednarz & Ostrowski (1998) also calculated 
the index adopting the small angle scattering model, which 
is approximately equivalent to the direction angle diffusion 
model in a mathematical sense. These works showed the 
spectral index a typically becomes 2.2 for ultra-relativistic 
shock limit, while our large angle scattering model gives 
<y ~ 1. This index of 2.2 may be universal when the direc- 
tion angle diffusion (or equivalent process) is proper for the 
scattering/deflection process in both upstream and down- 
stream, though it is not yet settled whether the real scatter- 
ing/deflection process during such quite short time can be 
treated well as a diffusion in direction angle. 

Because our probability function based description of 
the shock acceleration described in Section 5.1 needs only 
two probability functions Puw and Pdw to obtain the spec- 
tral index, it can be extended to include other scattering 
or deflection mechanisms if corresponding probability func- 
tions can be defined. For example, for the model of regular 
deflection in upstream the probability function Puw can be 
denned and expressed in an analytical form under the con- 
dition that the direction of velocity perpendicular to the 
shock normal is randomly distributed at the shock crossing 
from downstream. Thus, if the large angle scattering model 
is adopted in downstream (assuming the Bohm diffusion), 
one can calculate the spectral index. 



7 CONCLUSIONS 

We have given a new formulation of the first order Fermi 
acceleration in shock waves with any shock speed in the 
point of view of the random walk of single particles suffer- 
ing from large angle scattering. First, we have investigated 
the properties of particle trajectories in a moving medium. 
We showed that the problem could be formulated in terms 
of the random walk with absorbing boundaries based on 



the probability theory and derived an integral equation for 
the density of scattering points. By approximately solving 
it in an analytical form, we have given approximate ana- 
lytic expressions for the probability density of pitch angle at 
return for particles which cross the shock front at a given 
pitch angle. We have confirmed that our approximate solu- 
tions agree with the Monte Carlo results quite well for the 
isotropic scattering. We have also shown that they are quite 
different from those based on the diffusion approximation 
for relativistic and mildly relativistic shocks. It is seen that 
the non-diffusive effects such as the return after only a few 
steps of scattering or the finite mean free path, which are not 
included in the diffusion approximation, are very important. 

Using these results we have applied our formulation to 
the shock acceleration and compared our results with those 
in the literature such as 'the relativistic diffus ion approxima- 
tion', which was developed by Peacock (1981), together with 
the conventional diffusion approximation. We have found 
that the multi-step approximation, in which the effects of 
the return after only a few steps of scattering are neglected, 
gives a quite good fit to the Monte Carlo results on the pitch 
angle distributions at the shock crossing and the spectral 
index of accelerated particles. This is somewhat surprising 
because the multi-step approximation gives a poor fit to the 
Monte Carlo results when we fix the initial pitch angle at the 
shock crossing. One reason for this is considered that when 
averaging over the initial pitch angle, the effects of a few 
steps of return and the finite mean free path tend to cancel 
out and as a result the diffusive return term becomes domi- 
nant. Thus, our approach gives a theoretical base to use the 
multi-step approximation in the shock acceleration with any 
shock speed. We have also given a correct expression for the 
diffusive length scale for relativistic shocks, whic h equ als to 
the scale lengths previously derived by Peacock ( 1981 ) and 
Kirk & Schneider (1988) for far upstream distribution in rel- 
ativistic shocks, instead of a naive diffusion length based on 
the conventional diffusion equation. 
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APPENDIX A: THE RELATION BETWEEN 
THE DENSITY OF SCATTERING POINTS AND 
THE PHYSICAL NUMBER DENSITY 

In order to consider the connection of the density of scatter- 
ing points in random walk to the physical number density, 
let us introduce the concept of the staying time density. We 
define the staying time density S(x)dx as the total staying 
time in the region between x and x+dx for the particle which 
is injected at the origin x = and absorbed eventually by 
one of the absorbing barriers. The total staying time in the 
non-absorbing region before absorption is clearly given by 



S(x)dx. 



Then, the particle-averaged staying time density S(x) (i.e., 
the ensemble average of S(x)) is given as follows. Consider 
the contribution to the staying time of a particle between the 
successive scatterings. Since the probability of the particles 
to move beyond Ax without scattering is exp( — \Ax\/\* (/j,')) 
and since the staying time density is the pitch angle 

averaged staying time density made by the particle during 
one step, T(Ax), is given by 



T(Ax) 




P , — e **W)dii {Ax > 0) 



P,\ — \e^TPl dfi' (Ax<0) 
i v t l 



(Al) 



where v^i is measured in the boundary rest frame. Using the 
probability density of the m-th scattering point f m (x), the 
staying time density made by the particles during the m-th 
step to m + 1-th step is written as 



T(x — x')f m (x')dx' 



Therefore, summing up all steps from m — to oo, the total 
staying time density is written as 

S(x) = f T(x - x')n(x')dx + T(x). 

J -b 

If we consider a steady state problem, in which particles are 
injected at the origin at a constant rate (with pitch angle 



distribution P M >), S(x) gives the physical number density for 
unit injection rate. These relations also have been used in 
the Monte Carlo simulations (e.g. Naito & Takahara 1995). 



APPENDIX B: DERIVATION OF AN 
ALTERNATIVE FORM OF THE INTEGRAL 
EQUATION 

The problem of the random walk of single particles, which is 
discussed in Section is equivalent to a steady state prob- 
lem in which one particle is injected at the origin per one 
step. In the latter problem, f m (x) means the number density 
of scattering points of particles injected before m-steps and 
n(x) correspond to the number density of scattering points 
made by all particles existing that time. Here, we consider 
the integral equation for n(x) in this equivalent steady prob- 
lem. 

Let us consider the random walk of particles in steady 
state with no boundary but with imaginary boundaries at 
x — a and x — —b. The number density of scattering points 
clearly equal to njv(x). Then, we can distinguish the parti- 
cles that cross the boundaries at least once ('passed parti- 
cles') from those that never crossed the boundaries ('non- 
passed particles'). We write densities of scattering points 
of the former and latter particles as n pass (a:) and n(x), re- 
spectively. This n(x) obviously coincides with the density 
of scattering points of the random walk of single particle 
with the absorbing barriers at x = a and x — —b. The 
sum of these two densities should be the solution for the 
no-boundary problem, un(x), 



n N {x) = n(x) + n paBB (x). 



(Bl) 



Non-passed particles are injected at the origin steadily at 
one particle per step. On the other hand, passed particles are 
created when non-passed particles cross one of the imaginary 
boundaries and this can be regarded as the passed particles 
are injected at the first scattering point after they first cross 
the boundary. If we denote density of the first scattering 
point as D(X), this is calculated by the distribution of non- 
passed particles n(x) as, 



D(X) = f(X)+ / f(X-x )n(x )dx' (X > a or X < -6).(B2) 
J-b 

Each passed particle which is injected at each first scattering 
point Xi n j makes the density of scattering point un{x — 
Xi n j). Thus, n pass (x) is written as 



^pass (^e) 



n N (x~ X)D(X)dX 



n N (x- X)D{X)dX. 



(B3) 



Using equation (Bl) and (B3), we obtain 

f> oo 

n(x) = n N (x)- j n N (x - X)D(X)dX 



n N {x- X)D(X)dX. 



(B4) 
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By the definition of D(X), 

n N (x~ X)D{X)dX = / n N (x - X)f(X)dX 

J a 

pa poo 

+ / n(x') / n N (x-X)f(X-x')dXdx, (B5) 



-b p-b 

n N (x - X)D(X)dX = / n N (x- X)f(X)dX 



+ / n(x') / n N (x - X)f(X - x )dXdx'. (B6) 

J — b <j — oo 

Thus, defining the kernel as 

/■OC 

fe(x,a;') := / n N (x - X)f(X - x')dX 



+ / n N (x-X)f(X-x')dX, (B7) 

■/ X. 

we get a Fredholm integral equation of the second kind, 

r 

n(x) = un{x) — k(x, 0) — / k(x,x')n(x')dx' . (B8) 

This provides an alternative form of the integral equation 
for n(x). 



where we define functions 

V(b) := ^!{(A.(-^ + l + p) + /9 )e(6) 



+(— \- +A + p)R_E 2 (^-) + rLe 3 (^-) 
. , R- , „ , 6 



6 

Ir7 



:= E 1 (±)-e-^5E 1 (^- 
and a constant 



A' := no (ft- + g+(L )L ) + Ail 



"o 



2 1 -R- 



(D2) 
(D3) 

(D4) 



APPENDIX C: DERIVATION OF INTEGRAL 
CONDITION FOR N+(X;p,' ) 

This is the case of a — » oo. We first rewrite the Wald's 
identity ( |l9[ ) as 
6 



/(X - x) [n{x'; oo, 6) + 5(a/)} dx'dX = 1. 



This equation can be rewritten further as 



f(X - as') {n(x - 6; oo, 6) + 5{x' - b)} dx'dX 



- e ° b .(Cl) 



On the other hand, N+(x; h'q) in the equation ( p9| ) can be 
rewritten as 

1 f°° -At 
N+(x;fj,' ) = — / {n{x -b;oo,b) + S(x -b)}e x o db. 



o Jo 

Thus, combining these two, we obtain the condition to be 
satisfied for N+(x; n'o): 

f(X-x')N+(x';iJ,' )dx'dX = 1 (C2) 
1 + y oA 

Using the representation of p.d.f. for large angle scattering 
model (7), this condition becomes equation (143) . 



APPENDIX D: THE ABSORPTION 
PROBABILITY BY OUR APPROXIMATION 

The absorption probability from initial position 6, P a bs(&), 
by our approximation is calculated as follows. 



Paha{b) = A'e L o + v (b) 



(Dl) 
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